Quantum cooling by unitary transformations 



Winny O 'Kelly de Galway^ and Jan Naudts^ 

1 Instituut voor Theoretische Fysica, Universiteit Leuven, K.U. Leuven, B-3001 
! Leuven, Belgium 

2 Departement Fysica, Universiteit Antwerpen, Universiteitsplein 1, B-2610 



(N 



^ . Antwerpen, Belgium 

' E-mail: winny.okellydegalwayOfys.kuleuven.be, jan.naudts@ua.ac.be 

l> . 

Abstract. We study the unitary time evolution of a simple quantum Hamiltonian 
(— I ' describing a heat engine coupled to two heat baths. The engine is modeled as a three- 

Qh' level system. Each heat bath consists of a single harmonic oscillator. The engine 

(-H I is operated via time-dependent external fields. The 5-matrix of the thermodynamic 

cycle is obtained in analytic form. We conjecture that the spectrum of this S'-matrix 
contains a continuous part, and that this is a requirement for the operation as a heat 
engine. Energy currents flow in both directions through the engine. The balance of 
these currents determines whether the engine performs work or whether its operation 
^ , requires the application of external forces. 

in 
in 

(N 



1. Introduction 



Quantum heat engines produce work or pump energy by repeated execution of a 
thermodynamic cycle, much in analogy with the original Carnot engine or the Otto 
engine [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13]. The theoretical modelling of these engines 
^ I usually relies on the following ingredients. The state of the engine is described by a 
5^ I density matrix. During phases 1 and 3 of the cycle the engine is in contact with one of 
the heat baths. The time evolution of the density matrix is then described by a master 
equation of the Lindblad form [14, 15]. During phases 2 and 4 the state of the engine 
is modified by external forces working upon it. The adiabatic theorem is invoked to 
obtain that energy levels are modified without changing their occupational probability. 

More detailed models going beyond this standard approach are needed to under- 
stand some fundamental questions about topics such as the entanglement between the 
quantum engine and the heat baths or the role of processes of decoherence in the heat 
baths. An early attempt [16] in this direction, using the exact solution of a model 
describing a system in interaction with a heat bath, reveals important deviations when 
the decoherence processes are failing. 

It is rather straightforward to improve the modelling of the heat engine during the 
operation of the engine (phases 2 and 4). One of the techniques of quantum cooling 
is known as the STIRAP technique [17, 18, 19]. It involves population inversion of 
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three level systems. Detailed theoretical studies are available on how the three level 
system behaves when two nearly simultaneous light pulses are applied [20, 21]. Also 
analytical results are available [22] describing the behaviour of three level systems in 
time-dependent external fields. 

It is less easy to improve the description of the interactions between the engine and 
the heat baths. A full quantum mechanical treatment of a large heat bath seems not 
to be feasible. Never the less some results have been obtained by exphcitly treating the 
heat bath [23, 24, 25]. We follow this line of research and limit each of the two heat 
baths to a single harmonic oscillator. 

An additional motivation for our study comes from the investigation of thermo- 
dynamics in small closed quantum systems [26, 23, 27]. It is obvious to say that one 
harmonic oscillator is warmer than the other if the initial conditions are such that the 
engine performs work while energy flows from warm to cold. This can then be used to 
define the temperature of each of the two heat baths. Our present results only form a 
first step of such a rather challenging programme. 

We use units in which h—1. 

The model is introduced in the next Section. The (S-matrix approach is explained 
in Section 3. The analytic expression for the (S'-matrix corresponding with one 
thermodynamic cycle is obtained. In Section 4 we analyse our results. Final conclusions 
and an outlook follow in Section 5. The details of our calculations are explained in the 
Appendices. 

2. The model 

The model Hamiltonian H consists of an unperturbed part Hq describing the cold and 
warm heat baths and the heat engine, to which are added time-dependent external fields 
operating the heat engine and time-dependent interactions between the cold and warm 
baths and the engine. 

It is quite common to model a heat bath with a collection of harmonic oscillators — 
see for instance [16] . To keep the model simple the cold and the warm heat bath are each 
made up by a single harmonic oscillator. All together, the unperturbed Hamiltonian 
reads 

Ho = ujia^a + Hgef + oj^c'c (1) 

The operators a and c are the annihilation operators of the cold and of the warm 
harmonic oscillator, respectively. The hamiltonian of the heat engine H2 is discussed 
below. 

The use of a three level system as a heat engine has been proposed half a century 
ago by Scovil and Schulz-Dubois [1, 2]. In particular, the Stimulated Raman Adiabatic 
Passage (STIRAP) technique [18, 20, 21] has become a very efficient tool [19] to change 
the level populations of molecular systems. In the present study we limit our selves to 
a rather primitive sequence of two square pulses, although more realistic pulses can be 
treated analytically as well [22]. 
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Hgei = 



v 



(2) 



The Hamiltonian of the three level system is given by 

/ \ 

IX . 

n + J 

The three levels are labelled g, e, and /, and have energies —/i, /i, and ii+25, respectively. 
To this three level system we add two time-dependent external fields which are used to 
pump the energy levels of the heat engine. Their contribution is 

-^gef = -ea{t)Xl - eb{t)X6, (3) 

where Ai and Xq are the Gell-Mann matrices — see the Appendix A. 

The interaction between the three level system and each of the harmonic oscillators 
is inspired by the Jaynes-Cummings model. The choice for this model is in the first place 
justified by the consideration that its eigenvalues and eigenvectors can be calculated 
analytically. 

The coupling at the cold side is given byf 

Hi2^ Ki2{t){a^E+ + aE_) 



(4) 



with 



E+ = -(Ai + iAa) 




and 



(5) 



It couples the g and e levels of the three level system. At the warm side the interaction 
Hamiltonian is given by 

H23 ^ K23{t){F+C^ + F^C) 



(6) 



with 



/ \ 



^^(Xe + iXj) 



■{Xq - iXr) 



1 



/ \ 



1 



and 







(7) 



It couples the e and / levels of the three level system with the hot harmonic oscillator. 
The total time-dependent Hamiltonian is now 

H^Ho + H,2 + Isei + H23. (8) 

I In the Jayncs-Cummings model is multiplied with cr+ instead of C7_. The change made here is 
needed because the ground state of our three level system corresponds with the excited state in the 
Jaynes-Cummings model. 
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3. Cycles 

The external field strengths ea(t) and eb(t) and the coupling parameters nuit) and 
K23it) all depend on time t. They are pulsed one after another in such a way that a 
thermodynamic cycle is traversed. See the Figure 1. 




Figure 1. The 4 phases of the thermodynamic cycle. On the horizontal axis is the 
energy of the heat engine. On the vertical axis is the total energy of the system. 

The cycle starts by coupling the heat engine to the cold bath. The switching on and 
off changes the total energy of the system (this contribution is omitted in the figure). 
But during the first phase of the cycle the total energy is constant. In the second phase 
the energy of the heat engine is pumped up by applying a sequence of two pulses. Work 
is performed by doing so. In phase 3 the heat engine releases energy to the warm bath. 
In phase 4 the engine delivers work to the environment. This is again modelled by two 
externally applied pulses which pump down the internal energy of the heat engine. 

Note that the cycle does not necessarily close. It is obvious that in the energy 
transfer mode the engine will consume more work than it can deliver. Because the heat 
baths are finite this implies that the total energy of the system goes up after every cycle 
of the process. The time inversion symmetry of the purely mechanical system is broken 
by the sequence in which external pulses are given and couplings are activated. See the 
Figure 2. 
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Figure 2. Overview of the activation of the time-dependent terms in the Hamiltonian. 



Quantum cooling by unitary transformations 



5 



The S-matrix approach 

The time evolution of the system with Hamiltonian (8) is studied without making any 
approximation. This is best done in the interaction picture. Then the wave function 
of the total system — heat engine plus baths — is time-independent in the periods 
when none of the time- dependent terms is active. The effect of activating one of the 
interaction terms or one of the external fields is then to transform the wave function ijj 
by means of an S-matrix into a new wave function Sip. 



Step 1: Absorbing energy from the cold system 



In the first phase of the cycle the three level system is connected to the cold harmonic 
oscillator during a time ri. The corresponding S-matrix is denoted 5*1. It is not very 
difficult to calculate it exactly. See the Appendix B. The result is of the form 



with 



and 



Si 



^iriHo^—iTiH 

gin(^i-2M) ^^t(^ _ iC)aEi + ia.^BE+] 
+ e" 3-1 (-'1-2/^) [iBaE_ + aa^A + iCj^s] 
-I- GiEi + E^, 



(9) 
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(10) 



A 
B 
C 



V - cos (ti A„) I n) (n I , 
^-^ n-\-l 

n 

yZ I T- sin(ri A„) sin(26'„) \n) (n| , 

^-^ y/n -\- 1 

n * 

— — sin(riA„) cos(26'„)|n) (n|. 
^-^ n-\-l 



The coefficients A„ and the angles dn are given by 

1 



An 

tan(6ln) 



(12) 



2An + cui — 2^ 

The operator G\ is the orthogonal projection |0)(0| onto the ground state of the cold 
harmonic oscillator. 



Quantum cooling by unitary transformations 



6 



Step 2: Pumping up 



We apply a sequence of two pulses of the on/off type. The first pulse results by giving 
ea{t) a constant non-zero value during a time Ta- It tries to invert the population of the 
levels e and /. The change of the population as a consequence of this pulse is given by 
the S-matrix which is now calculated. 



jTaH, 



(13) 



Note that we switched notations, using two-dimensional Pauli matrices instead of the 
Gell-Mann matrices, omitting one dimension for a moment. Introduce the constant 
Ta = There follows 

S2a = [cOs{TaS) - ism{TaS)a3\ 

X [cos(t„/7;) + iTaSm{Ta/Ta){Sa3 - e„(7i)] . 

(14) 

Let us now make an appropriate choice of the pulse duration Ta- The goal is to 
minimize the population of the e-level after the pulse. Since one can expect that before 
the pulse the e-level is more populated than the /-level the best one can do is to require 
that the e matrix element of S2a is as small as possible in modulus. Let therefore 
Ta — \T^Ta. Then the S-matrix becomes 

S2a = iTa [cOs(t„(5) - 1 Sin(Ta(5)(73] SUl{Ta/Ta){5a2, - ea(7i). (15) 

Restoring the third dimension this becomes 
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+ iTa 




5e-'^-^ 
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'J 
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-Sae'^-^ 





In the limit of a strong short pulse this becomes 

/ 1 



■>2a 



\ 




-i 



(16) 



(17) 



The first pulse of the second phase of the thermodynamic cycle is followed by a pulse 
of duration r^, intended to invert the population of levels e and g. The corresponding 
S-matrix reads, using the notation Tj, — I/a/a*^ + ef, 



2b 



= [cos(Tb//) - isin(Tb/i)(73] 

X [cos(T6/r6) + iTbSm{Tb/Tb){iia3 + ebai] 



(18) 
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With similar arguments as before let us choose tj, — hnTi,. Then the S-matrix becomes 





/ 















/ 













'S'26 — 















+ iTb 













(19) 
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In the limit of a strong short pulse this becomes 

/ i \ 



'26 



v 



i 
1 



(20) 



All together the S-matrix for the second phase of the cycle equals 

iTbtbe'^'''' 




S2bS2a 

( iTh, 



\ 



iTaS 



-TaTbiieae-'-'^^+'-»^' 



In the limit of strong short pulses it becomes 
/ 1 \ 



i 
-i 



(21) 



(22) 



Step 3: Exchanging energy with the warm bath 

In the third phase of the cycle the three level system is connected to the warm harmonic 
oscillator during a time T3. The corresponding S-matrix is denoted S3. The calculation 
is similar to that in Step 1. The result is of the form 



S. 



gi73-ff0g-2T3-ff 



— El + E2G3 

+ e^^3(<^3-25) ^gt(^ _ iV)cE2 + ic^YF+] 
+ 6-^^3(0.3-25) ^YcF_ + cc\Z + iV)E3] 



(23) 



with 



1 



cos(T3^„)|n)(n| 
sin(r3Cn) sin(20„) I n)(n I, 



n 

y = ^ — — - sin(T3^n) cos(20„) \n) {n\ . 
n I i. 

n 

The coefficients ^„ and the angles 0„ are given by 



(24) 
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tan(0„) 



(25) 



The operator is the orthogonal projection |0)(0| onto the ground state of the warm 
harmonic oscillator. 

Step 4- Pumping down 

The operation in the fourth phase is the inverse of that in the second phase. We thus 



4. Analysis 

In the previous Section the contribution to the S-matrix from each of the four phases 
of the cycle has been obtained. The composite matrix S — S4SSS2S1 is now calculated. 
The result is a rather complicated. Therefore a tensor notation is appropriate. Remem- 
ber that the Hilbert space of wave functions of the total system is the tensor product 



The first and the last factor are the Hilbert space of the cold and of the warm harmonic 
oscillator, respectively. The middle factor is the space of vectors with three complex 
components. 

4-1- The composed S-matrix 
The full S-matrix reads 



have S4 = Sl- 



•warm- 



(26) 



S^ 



X 



X 



X 



X 




(27) 
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Table 1. Interpretation of the terms appearing in (28). 
The arrows indicate the direction of the energy flow. 

at {A -iC)a® El ® c){Z-iV)c — — 



a'^B® E+®c'^{Z -iV)c ^ — 

Ba®Ei® c^Y -J> 

aa<{A-\-iC)®E+®c<Y — 

a){A-iC)a®E_®Yc — ^ 

a^B ®E2®Yc ^ <s- 

Ba®E_®cc\Z + iV) — 



aa<{A + iC)®E2®cc<{Z + iV) — — 



For simplicity, we use the value (22) of 5*2 in the limit of strong short pulses. In this limit 
one has ^2^^152 = ^3, 5|^2^2 = ^i, SIe^S2 = ^2, SIF+S2 = -E+, SIf_S2 = -E_. 
Hence, the above expression for S simplifies to 

+ ei^3(<^3-25) ^i^Ei<^ c\Z - iV)c -iI^E+<^ c^Y] 

+ e-i^('^3-25) ^E^0Yc + I®E2® cc\Z + iV)] } 

X |e^i('^i-2/^) [a\A - iC)a ®Ei®I + ia^B ®E+®t\ 
+ e- ('^1-2'^) [iBa (g)E^0l + aa\A + iC) 0^2®!] 
+ Gi ® ^1 1 + 1 ® 1| 

= I ® ® I + Gi ® El ® G3 

+ e^^i^'^i-^") [a\A - iC)a ® ® G3 + ia^B ^ E+ ® G3] 

+ e^^^'^^-^^^Gi ® c\Z - iV)c - ie-5^3('^3-25)^^ ®E_®Yc 

_|_ gt75('<J3-2i5)g|ri(a)i-2/i) 

X [a\A - iC)a ®Ei + ia^B (g) E+] c\Z - iV)c 
^ gi.g(a,3-25)g-|n(^i-2M) ^Ei- iaa^{A + iC) ® E+] ® c^Y 
+ g-ir3(c.3-25)gfnK-2M) \_ia\A - iC)a ®E^ + a^B (g) E2] ® Yc 

_^ g-iT3(a;3-25)g-|Ti(a;i-2/.) 

X [iSa (g) -B_ + aa^{A + iC) ® E2] ® cc\Z + iV). 

(28) 

Note that the operators A, G, y, Z, l^, commute with the counting operators of the 
two harmonic oscillators. Hence the two terms which directly transfer energy between 
the two oscillators are those proportional to Ba®Ei®c^Y and a'^ B®E2®Y c respectively. 
They act in opposite directions. Other terms do not transfer energy or they exchange 
energy between the heat engine and one of the oscillators. See the Table 1. 

4-2. Eigenvectors of the S-matrix 

The above S-matrix describes the effect in the interaction picture of performing one 
thermodynamic cycle. It is immediately clear that the ground state |0, g, 0) of the system 
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is an eigenstate of this S-matrix with eigenvalue 1. This is an immediate consequence 
of the fact that the ground state of the Jayncs-Cummings model is not affected by the 
interactions of the model. An important question is whether the S-matrix has other 
eigenvectors. Indeed, such eigenvectors describe situations in which the action of the 
heat engine has no effect at all. Of course, on a superposition of eigenvectors the engine 
can have effect. But the result is an almost periodic function which always returns 
arbitrary close to its starting point. 

On the other hand, if part of the spectrum of S is continuous then a genuine energy 
transfer is possible by which the system approaches a stationary regime. We did not 
yet succeed to prove the existence of such a continuous spectrum. But note that several 
terms in the expression (28) are linear in the creation and annihilation operators of 
the cold and warm oscillators. These operators are linear combinations of position and 
momentum operators which are known to have a purely continuous spectrum. Hence it 
would not be strange when this continuous spectrum survives when adding other terms. 



4-3. Energy transfer 

The result (28) seems hopelessly complicated but can never the less be used to derive 
some unexpected properties of the heat engine. The change in the state of the cold 
oscillator before and after one cycle is defined by 

D = S^a^aS - a^a. (29) 
One finds (see the Appendix C) 

D = aa^B'^ ^ E2 - a^'B'^a ® Ei 

+ ia^aa^{A + tC)B ® E+ -i{A-tC)Baa^a® E_ ® I. (30) 
The eigenvectors of D are of the form 

ijj ^ u\n + 1) (S> \g) + v\n) <^ \e) (31) 

(we neglect the Hilbert space of the warm oscillator for a moment). The condition 
Dip — fnp then yields 

0^u{p + sin2(riA„) sin^(26'„)) 

+ ■i;sin(TiA„) sin(26'n) [sin(TiA„) cos(26'„) - icos(TiA„)] 
= Msin(TiA„) sin(2^„) [sin(riA„) cos(2^„) + icos(TiA„)] 

+ v{p- sin'(TiA„) sm\20„)) . (32) 

This set of equations has a non-trivial solution when 

p = ±sin(riA„) sin(26'„). (33) 

Corresponding eigenvectors are then 

u — sin(riAn) cos(2^^„) — i cos(tiA„), 

^1- sin(riA„) sin(2^„). (34) 
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Note that -D|0) (g) \g) = 0. Hence, the spectrum of D is completely known. For 
each strictly positive eigenvalue p > also — p is an eigenvalue, p > corresponds with 
warming up of the left oscillator, p < with cooling. 

One concludes that raising or lowering the energy of the cold oscillator after one 
cycle of the engine depends completely on the choice of the initial wave function. The 
important question is of course what happens after one cycle with a wave function 
originally chosen as an eigenvector t/j oi D with negative eigenvalue. Will St/j be a 
superposition of eigenvectors all with negative eigenvalues? Or will part of them have a 
positive eigenvalue? We expect, but did not succeed to prove, that the latter is the case 
and that in this way the cooling efficiency of the engine gradually goes down so that 
after a small number of cycles a stationary regime is approached. 

A similar calculation for the warm oscillator is possible. But note that an easy 
result only follows when starting the cycle with coupling the heat bath to the warm 
oscillator instead of the cold one, as used in the above calculations. 

4.4- Performing work 

The previous subsection gives a partial answer to the question whether the engine is 
capable of transferring energy between the two oscillators. Now follows a discussion of 
the work needed to operate the engine. 

In phases 1 and 3 of the cycle some work is needed to operate the valves connecting 
the engine with the cold respectively the warm oscillator. Indeed, switching on and 
off the interaction terms (4, 6) changes the total energy of the system. Since the wave 
function of the system evolves in time between the switching on and switching off the 
involved energy changes to not necessarily cancel. Hence we expect that a tiny amount 
of work is needed to operate these valves. 

It is now indicated to consider a cycle starting with phase 2 instead of phase 1. 
Then the energy changes during the respective phases are given by 



AE, = Ho- SiHoSl 
AE2 = SIH0S2 — Hq, 
AEs = SliSlHoSs-Ho)S2, 
AEi = Slsl{S2HoSl - Ho)S3S2. 



(35) 



Using the simplified expression (17) for S2 one obtains 



AEi = (wi - 2p) {B^aa^E2 - a^B^aEi) 




(36) 



and 



(37) 



and 



AE3 = (^3 - 2(5) {E2CC^Y^ - Eic^Y^c] 
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- i{ujs - 28)E+cHZ + iV)cc^Y 
+ i{uj^ - 25)E^{Z - iV)cc^Yc 



(38) 



and 



d.Ei = - A£;2 + 2(// - 5) {E^c^Y'^c - E2Y^cc^) 

+ 2i(/x - S) [E+c\Z + iV)cc^Y - E_Ycc\Z - iV)c] , 



(39) 



where [a, 6, c] denotes the diagonal matrix with eigenvalues a, 6, c. See the Appendix D. 

Several features can be observed. The contributions A£'3 and A£^i represent the 
energy needed to switch on and off the interactions with the harmonic oscillators. They 
vanish when the coupling between the heat engine and the oscillators is at resonance. 

The work performed by the engine equals the quantum expectation of the operator 
—AE2 — AE4^. When fi = 6 then the operation of the engine is meaningless and no net 
energy is used and no net work is performed during the phases 2 and 4. In the general 
case the eigen values of AE2 + AE4 can be calculated analytically. One obtains 



The corresponding eigen vectors are linear combinations of \g, n+1) and |e, n) (neglecting 
the state of the cold oscillator). Hence also the spectrum of this operator is symmetric 
under a change of sign. This means that the initial conditions determine whether 
operating the engine consumes energy or whether it performs work. 

5. Conclusions 

It is feasible to obtain analytic results for a closed quantum system consisting of a heat 
engine operating between two small heat baths. The engine is operated by switching 
external fields on and off. The state of the system is at any moment determined by 
its wave function. The time evolution follows from solving the Schrodinger equation 
using a time-dependent Hamiltonian. In the traditional approach the engine is the 
system and the cold and warm heat baths belong to the environment which is taken 
into account in a phenomenological way. In the present approach the boundary between 
the system and the environment is shifted so that the heat baths belong to the system 
while the operation of the engine is still implemented in a phenomenological way. The 
main advantages of this approach are that the time evolution is unitary and that the 
quantum entanglement between the engine and the heat baths is treated rigorously. 
From our toy model we have learned a number of points. 

• The use of the interaction picture for doing the calculations improves the 
transparency of the calculations. 

• We do not make use of the adiabatic theorem. The change in the population of the 
energy levels of the engine results from the time evolution. As a consequence all 
results depend only on intra-level distances and not on the positioning of oscillator 
levels w.r.t. levels of the engine. 




(40) 
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• At each of the two interfaces the energy flows in both directions. Energy leaks away 
in the direction opposite to the intended one. Eight different energy contributions 
have been distinguished in the Table 1. In the usual approach these are replaced 
by two phenomenological terms. 

• We conjecture that the S'-matrix of the thermodynamic cycle has a continuous 
spectrum and that this is an essential feature of a realistic heat engine. The 
argument for the latter is that eigenstates of the ^'-matrix are unchanged after 
one cycle. Hence they correspond with states of the system in which the engine is 
inactive. 

• To our surprise the operator D — S^a^aS—a^a which measures the change in energy 
of the cold bath during one cycle has a fully discrete spectrum with explicitly known 
eigenvectors and eigenvalues. This is probably an artefact caused by the Jaynes- 
Cummings mechanism used by the toy model for implementing the interaction 
between heat bath and engine. However, the spectrum of this operator D is 
symmetric under the change of sign. This could be a more general feature being a 
consequence of time inversion symmetry. 

• The change of energy during one cycle can be obtained analytically as well. The 
operation of the valves connecting the engine with the heat baths costs energy 
except when the interaction is at resonance. The pumping up and down of the 
occupational probabilities of the engine levels can cost energy or can perform work 
depending on the initial state of the system, this is, depending on its wave function. 
This shows that the engine can be used either to transfer energy from the cold to 
the warm oscillator or to perform work produced by the energy flow from warm to 
cold. 

The present paper shows that our approach of studying a heat engine by means of 
unitary time evolutions can produce new insights. Future work goes into two directions. 

Many questions about the model are still unresolved. We have studied only one 
cycle of the engine. The repeated operation of the engine as an energy pump will lead 
to a steady increase of the total energy. But the reduced density operator of the cold 
oscillator might converge. Alternatively one could try to analyse how the effect of a 
single cycle depends on the initial wave function. 

It is tempting to analyse our results in terms of thermodynamic concepts. This is 
not easy because in a closed system it is not clear at all what the distinction is between 
energy and heat. But one could for instance analyse the change of entropy of the reduced 
density matrices. 

Another line of research could try to study some features of the toy model in a 
more general context. Particularly intriguing is the question whether the continuous 
spectrum of the S'-matrix is a general feature of heat engines. Also the effects of the 
time inversion symmetry of the equations of motion require further investigation. 
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Appendix A. The Gell-Mann matrices 



Conventionally, the Gell-Mann matrices are defined as follows. 



Ai 



A, 



A, 




A, 



Ap 
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/o 



1 

/o 






As - ^ 
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/l 





-i 0\ 



1 \ 
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2/ 



(A.l) 



Appendix B. The S-matrix of phase 1 of the cycle 

Here we calculate the S-matrix of a Jaynes-Cummings Hamiltonian in which the 
interaction is switched on during a finite time. The coupling is constant with strength 
Ki2 during a time interval of length ti. The relevant Hamiltonian is 

H = ooia^a — iJ,az — K.uici'cr^ + aa-). (B-l) 

Let 



\9) 



1 




and |e) 




1 



(B.2) 



The eigenstatcs of the harmonic oscillator are denoted n), with n — 0,1, ■ ■ -. 
The ground state of H is 

|0,-) = |0)®b). (B.3) 

It satisfies H\0,—) — — /i|0, — ). The pairs of excited states are denoted \n,±), with 
n = 0, 1, • • •. They are of the form 



\n, -) = cos(6'n)|n) (g) |e) + sin(6'n)|n -hi)® \g), 
\n,+) = - sin(6l„)|n) (g) |e) + cos{9n)\n + 1) (8) \g). 



Prom 



follows 



H\n) (g) |e) = {ntoi + ii)\n) (g) |e) - Kv2,\/n + \\n 1) (g \g) 
H\n) (g \g) = (ncji - iJ,)\n) (g l^f) - KuV^ln - 1) g) |e) 

if |n, -) = cos(6'n) {ncoi + jj) \n) (g |e) - Ki2\/n + l\n -M) (g 



(B.4) 



(B.5) 
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+ sm(^„) {{n + l)uji - //) |n + 1) (g) \g) — K,i2Vn + l\n) 
H\n, +) = - sin(6'„) {nooi + /i) \n) (8 |e) - Ki2Vn + l\n + 1) (g) | 
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+ cos(^„) ((n + l)uji - /x) |n + 1) \g) - Ki2Vn + l\n) ® |e) 



The requirement that H\n,±) — E^\n, ±) then yields the set of equations 

{nui + /i) cos(^„) - Ki2Vn + l sin(6l„) = E~ cos(6'„) 

-Ki2Vn+ lcos(6'.„) + ((n + l)uji - n) sin(6'„) = 



sm 



= -E+ 



-Hi2Vn+ lcos(6'„) - {nui + fi) sin^ 
((n + l)uji - jji) cos{9n) + Ki2\/n + lsin(6'„) = cos(6'r, 



The solution is 



tan^„ 



n 



2A„ + 2// - a;i 
2Ki2\/n + 1 



with 



'^12(^ + 1) + ( A* - 0*^1 



A short calculation now gives 

Si\n) O |e) 



e^*^"e"**^|n) ® |e) 



COS 



)e-**^"|n,-) 



sm 



)e-**^" |n, +) 



cos(tA„) + icos(26'„) sin(tA„)] \n) |e) 



+ ie-^*("/2) sin(2^„) sin(a„)|n + 1) ® |^) 
+ 1) 1^) = e^*^°e-^*'^|n + 1) |^) 



= e 



COS 



|n, +) + sin(^„)e-**^" |n, -) 
•git(M-u,i/2) sin(2^„) sin(tA„)|n) ® |e) 
+ e-^*("/2) [cos(a„) - icos(2^„) sin(tA„)] 
X |n + 1) (g) 1^). 



(B.6) 



(B.7) 



(B.8) 



(B.9) 



(B.IO) 



These expressions can be written as (9). 



Appendix C. Change in the state of the cold oscillator 

Here we calculate (30). 

Note that [a'' a, Gi] = and [a'^a, a'^Aa] = and [a^a, a'^B] = a'^B and [a'^a, a'^Co] 
0. Using these relations one obtains 

[ata, S] = ie^^'^'^'-^f'^a^B ® ® G3 

+ ^e^3(^3-25)giri(u,i-2/.)^t5 ®E+(g) c\Z - iV)c 
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_ gir3(<^3-25)g-lriK-2M)^^ ® Ei ® C^Y 
^ g-iT3(a;3-25)giri(^l-2A*)^t5 ® ^2 ® 

- ie-^^3(a.3-25)g-fn(a;i-2M)^^ ^ E_ ^ Cc\Z + iV). (C.l) 

Note that 

X = ccHz'' + V^)+Y'' ^—^\n){n\ (C.2) 

n + 1 

n 

SO that cc^X — + c^Xc = I. Hence, from (C.l) one obtains 

D^S^[a^a, S] 

= + [ia.^aaHA + iC)B®E+ + aa^ B"^ ® E2] ® (Gg + c^Xc) 
- [a^B'^a ®Ei+ i{A - iC)aa^ Ba (g) EJ\ ® cc^ X 

= aa^B"^ ®E2- a^B'^a (g) Ei 

+ m^aa^(A + iC)B E+ - i{A - iC)Baa^a ® (g) I. (C.3) 

This is (30). 

Appendix D. Work performed during phases 1, 3, 4 

We first calculate AEi = Hq — SiHqSI. Note that one can write AEi — Si[SI,Hq]. 
Therefore we start with calculating 

[Sl Ho] = - i{uji - 2/i)e-^^(^^-2^)5a ® E_ 

+ ilui - 2/i)e^"i(^^-2^)a^5 ® E+. (D.l) 

Now multiplying from the left with 5*1 yields (36). 
Next calculate A£'3 using 

AEs = sIHoSs -Ho^ Sl[Ho, -S3]. (D.2) 

One calculates using [c^c, c\Z - iV)c] = [c^c, cc\Z + iV)] = [c^c, = 

[//o, = 5[E^ - E2, Ss] + u^ic^c, 83] 

= i{uj3 - 2S)e^^^'^'-^^^F+c^Y - i{uj3 - 2S)e-^^^'^'-^^^ F_Yc. (D.3) 

This gives using F_F+ = E3, E2F+ = F+, F+F_ = E2, E^F. = F_ 

SI[Hq, S^] = {003 - 26)EsYcc^Y + ^(^3 - 26)F+c^{Z + iV)cc^Y 

- {uj3 - 2S)E2c'Y^c - i{uj3 - 26)F_{Z - iV)cc^Yc. (D.4) 

It is then straightforward to obtain 

SI{SIHqS3 - Hq)S2 = (cus - 26)E2Ycc^Y - {us - 26)Eic^Y^c 

-i{ujs- 26)E+c\Z + iV)cc^Y 
+ i{uj3 - 2S)E_{Z - iV)cc^Yc. 

(D.5) 

This yields AE3. 

Finally calculate A£^4. One has using the simplified expression (17) 

S2HoSl-Ho = 2[ii + 6,-ii,-5]. (D.6) 
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Note that (using E2S2 = iE_, F+S2 = -iE2, and E3S2 = -iF_) 

S3S2 = I (8) E1S2 <^l + iI<S> E_(S>G3 

+ e5^3('^3-2^)l ® [iE_ ^c\Z- iV)c + E2® c^Y] 

+ e-i^3(o.3-25)j ^ ®Yc- iF_ ® cc\Z + iV)] . 

This gives (using SIE1S2 = £^3 and 4^2-52 = £^1) 

S2S\EiS3S2 = E3 



(D.7) 



SISIE2S3S2 ^I^ Ei^ Gs + I^ Ei^c^Z^ + V^)cc^c + I^E2^ Y^cc' 

-il(^ E+(^c^{Z + iV)cc^Y 

+ iI(»E_^ Ycc\Z - iV)c 
SISIE3S3S2 ^l®Ei® c^Y^c + 1®E2®{Z^ + V''){cc^f 

+ il® E+®c\Z + iV)cc^Y 

-il®E_® Ycc^Z - iV)c. (D.8) 



The result is 



AE4 = 2(11 + 5)E3 

- 2/iI ® [^1 ® G3 + ^1 ® c\Z^ + V^)cc^c + E2® Y^cc^] 

- 251 ® [E^ ® c^Y'^c + E2® {Z^ + V^){cc^f] 

+ 2i{^l - 5)1 ® E+® c\Z + iV)cc^Y -E_® Ycc\Z - iV)c . 

(D.9) 

Using G3 + c^Xc — 1 and cc^X — I this can be written as 

/\Ei = - A2 + 2(/i - 5)1 ®Ei® c^Y'^c - 2(/x - 5)1 (^E2<^ y^cc' 

+ 2i{fi - 5)1 ® E+® c\Z + iV)cc^Y -E_® Ycc\Z - iV)c . 

(D.IO) 

This is (39). 



References 

[1] H. E. D. Scovil and E. O. Schulz-DuBois, Three-Level Masers as Heat Engines, Phys. Rev. Lett. 
2, 262-263 (1959). 

[2] J. E. Geusic, E. O. Schulz-DuBois, and H. E. D. Scovil, Quantum Equivalent of the Carnot Cycle, 

Phys. Rev. 156, 343-351 (1967). 
[3] CM. Bender, D.C. Brody and B.K. Meister, Quantum mechanical Carnot engine, J. Phys. A33, 

4427-4436 (2000). 

[4] M. O. Scully, M. S. Zubairy, G.S. Agarwal, and H. Walthcr, Extracting Work from a Single Heat 
Bath via Vanishing Quantum Coherence, Science 299, 862-864 (2003). 

[5] T. D. Kieu, The Second Law, Maxwell's Demon, and Work Derivable from Quantum Heat Engines, 
Phys. Rev. Lett. 93, 140403 (2004). 

[6] T. D. Kieu, Quantum heat engines, the second law and Maxwell's daemon, Eur. Phys. J. D39, 
115-128 (2006). 

[7] Y. Rezek and R. Kosloff, Irreversible performance of a quantum harmonic heat engine. New J. 
Phys. 8, 83 (2006). 

[8] H. T. Quan, Yu-xi Liu, C. P. Sun, and F. Nori, Quantum thermodynamic cycles and quantum 
heat engines, Phys. Rev. E 76, 031105 (2007) 



Quantum cooling by unitary transformations 



18 



[9] H. T. Quan, Quantum thermodynamic cycles and quantum heat engines. II, Phys. Rev. E79, 
041129 (2009). 

[10] P. Salamon, K. H. HofFmann, Y. Rezek and R. KoslofF, Maximum work in minimum time from a 

conservative quantum system, Phys. Chem. Chem. Phys. 11, 1027-1032 (2009). 
[11] M. Esposito, R. Kawai, K. Lindcnbcrg, C. Van den Broeck, EfEciency at Maximum Power of 

Low-Dissipation Carnot Engines, Phys. Rev. Lett. 105, 150603 (2010). 
[12] M. Esposito, R. Kawai, K. Lindenberg, C. Van den Broeck, Quantum-dot Carnot engine at 

maximum power, Phys. Rev. E81, 041106 (2010). 
[13] J. Naudts and W. O'Kelly de Galway, On the BCH formula of Rezek and Kosloff, Physica A390, 

3317-3319 (2011). 

[14] A. Kossakowski, On quantum statistical mechanics of non-Hamiltonian systems, Rep. Math. Phys. 
3, 247-274 (1972). 

[15] G. Lindblad, On the generators of quantum dynamical semigroups, Commun. Math. Phys. 48, 
119-130 (1976). 

[16] N. G. van Kampen, A Soluble Model for Quantum Mechanical Dissipation, J. Stat. Phys. 78, 

299 310 (1995). 

[17] K. Winkler, F. Lang, G. Thalhammer, P. van der Straten, R. Grimm, J. H. Denschlag, Coherent 
Optical Transfer of Feshbach Molecules to a Lower Vibrational State, Phys. Rev. Lett. 98, 
043201 (2007). 

[18] K. Bergmann, H. Theuer, B. W. Shore, Coherent population transfer among quantum states of 

atoms and molecules. Rev. Mod. Phys. 70, 1003-1025 (1998). 
[19] E. Kuznetsova, M. Gacesa, Ph. Pellegrini, S. F. Yelin, R. Cote EfRcient formation of ground-state 

ultracold molecules via STIRAP from the continuum at a Feshbach resonance. New J. Phys. 11, 

055028 (2009). 

[20] V.N. Vitanov, Analytic model of a three-state system driven by two laser pulses on two-photon 

resonance, J. Phys. B31, 709-725 (1998). 
[21] K. Na, L.E. Reichl, Nonlinear dynamics of ladder and lambda STIRAP, Chaos, Sohtons & Fractals 

25, 185-196 (2005). 

[22] J. Naudts and W. O'Kelly de Galway, Analytic solutions for a three-level system in a time- 
dependent field, Physica D240, 542-545 (2011). 

[23] M. Campisi, P. Talkner, and P. Hanggi, Finite bath fluctuation theorem, Phys. Rev. E80, 031145 
(2009). 

[24] M. Campisi, P. Talkner, and P. Hanggi, Thermodynamics and fluctuation theorems for a strongly 
coupled open quantum system: an exactly solvahlc case. J. Phys. A42, 392002 (2009). 

[25] M. Esposito, K. Lindenberg, C. Van den Broeck, Entropy production as correlation between system 
and reservoir, New J. Phys. 12, 013013 (2010). 

[26] P. Talkner, P. Hanggi, and M. Morillo, Microcanonical quantum fluctuation theorems, Phys. Rev. 
E77, 051131 (2008). 

[27] M. Esposito, R. Kawai, K. Lindenberg, C. Van den Broeck, Finite-time thermodynamics for a 
single-level quantum dot, Eur. Phys. Lett. 89, 20003 (2010). 



